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Abstract 

Nanoscale objects often behave differently than their 'normal-sized' counterparts. Sometimes it is enough to be small 
in just one direction to exhibit unusual features. One example of such a phenomenon is a very specific in-plane 
magnetic anisotropy observed sometimes in very thin layers of various materials. Here we recall a peculiar form 
of the free energy functional nicely describing the experimental findings but completely irrelevant and thus never 
observed in larger objects. 
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1. Intriguing experimental observations 

In [1] we find the experimentally observed in-plane 
magnetic anisotropy energy (MAE) diagrams for mul- 
tilayer structure Cr(4)/Fe(2)/Cr(fl? Cr )/Fe(4)/Cr(2), where 
numbers are thicknesses of components, expressed 
in nm. The thickness of the middle Cr layer, dci, was 
varied in few nm range, and the complete structures 
were deposited on Si(100) substrate, covered with nat- 
ural SiC>2 layer 1. 5-2.0 nm thick. The substrate was 
not perfectly flat — as a result of ion beam erosion 
it was covered with quite well ordered ripples (see 
atomic force microscopy (AFM) images of the sub- 
strates, Fig. 1 in [1]). The metallic layers were deposited 
using molecular beam epitaxy (MBE) technique. Then- 
transmission electron microscopy (TEM) cross sections 
revealed mostly amorphous structure with small inclu- 
sions of polycrystalline character. The values of MAE 
were derived from hysteresis loop area observed while 
exciting field was oriented along successive in-plane di- 
rections. 

In principle, samples of this kind should not exhibit 
any in-plane magnetic anisotropy. This is indeed the 
case when the substrate is flat (see Fig. 4a in [1]). On the 
rippled substrate, however, this is no longer true and the 
sample exhibits peculiar two-fold in-plane anisotropy 
(coercive field, Fig. 4b in [1], MAE — in Fig. 4c). It 
is peculiar since it is not the uniaxial anisotropy: four 
maxima are visible instead of just two. 
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2. The surface magnetic anisotropy of a cylinder 

Consider the static configuration of individual spins 
located on a surface of a long (ideally: infinitely long), 
hollow ferromagnetic cylinder. In absence of any exter- 
nal field one may expect that individual spins may adopt 
one of the three stable configurations: 

• they all may be aligned with Coo symmetry axis 
of a cylinder. This is the lowest exchange energy 
configuration. 

• they all may be oriented perpendicularly to the 
above symmetry axis. Now the exchange energy 
is no longer at its global minimum. Nevertheless, 
such a configuration is stable since it realizes a lo- 
cal minimum of exchange energy. 

In the second case we may again distinguish two cases: 
either individual spins are aligned with local C2 sym- 
metry axis (there are infinitely many of them, each per- 
pendicular to Coo) thus pointing inwards or outwards of 
a cylinder and being perpendicular to the cylinder's sur- 
face, or they can be perpendicular to the local C2 axis 
(laying on the cylinder's circumference), making a ring- 
shaped configuration, and producing no net magnetiza- 
tion. It is easy to see that both those configurations are 
energetically equivalent, since the angle between any 
two neighboring spins is exactly the same. 

Here we concentrate only on this part of magneto- 
static energy which originates from Heisenberg-type ex- 
change interactions between spins. For a pair of nearest 
neighbor (nri) spins, say i and j, we have 



E = -2 J § ' i ■ § j = -2JS 2 cos (fij, 



(1) 
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where J is an exchange integral, tpij is the angle between 
spins ;' and j, and S = = js^-j. As the angle tptj 
between neighboring spins is small, then the following 
approximation is valid: 

cosp y « 1-^/2 = > ( 2 ) 

where 5 is the spacing between nn spins, and p is the 
cylinder radius. 

Full magnetostatic energy of a sample is of course the 
sum, running over all the pairs of nn, of expressions like 
(1). Anisotropy characterizes differences of free energy 
between various directions of an external field, so any 
constant terms are meaningless and may be dropped. In 
our case such a term is "1" in (2). After this is done, the 
exchange energy for a single nn pair of spins reads: 

S 2 6 2 

E*J— (3) 

P 2 

As the sum of expressions of type (1) is hard to treat 
analytically, we replace it with appropriate integral, i.e. 
we assume the continuous distribution of interacting 
spins but we do not approximate anything else. Partic- 
ularly, we do not make use of approximation (2). This 
way we have to integrate proper expression along the el- 
liptical path, being a trace of a cylinder's cross-section 
by a plane parallel to the external magnetic field. The 
final result for the surface part of the free magnetostatic 
energy density, E s , already presented some time ago in 
[2], reads: 

E s = K s \cosO\ (4) 

Here K s is the surface anisotropy constant, and 6 de- 
notes the angle between the direction of sample's mag- 
netization and easy direction Coo. As one might expect 
K s oc J/p 2 - in full accordance with simplified ap- 
proach, sketched in (3). 

A comment is in order in this place. Magnetocrys- 
talline anisotropy energy density is always expressed by 
even powers of cos 8 and is always a smooth function 
of the external field orientation. Here we have cos 9 in 
first power, and, additionally, the energy density is not 
a smooth function. 

3. Experimental confirmation 

The formula (4) has been first derived to in- 
terpret ferromagnetic resonance (FMR) spectrum of 
Co68Mn 7 SiioBi5 glass-coated amorphous single mi- 
crowire with diameter roughly equal to 16.5 /an. The 



spectrum, taken at fixed frequency, and containing more 
that one absorption line, could not be described (mod- 
eled) satisfactorily with conventional two- and fourth- 
order uniaxial anisotropics alone [3]. Unfortunately, 
even the inclusion of the surface anisotropy term (4) 
into the full expression for the free energy density didn't 
help much. This applies also to further experiments, per- 
formed on similar but thinner wires, down to the diame- 
ter of 6 pm. Some qualitative features of the spectra (e.g. 
broadening and distortion of absorption lines at special 
orientations), however, could be attributed to the pres- 
ence of a non-smooth surface anisotropy term. Never- 
theless, it had to be concluded that the wire's diame- 
ter was most likely to big to clearly observe the surface 
anisotropy contribution. By the way, due to the pres- 
ence of a glassy cover, other effects, notably the mag- 
netostriction of the inhomogeneously stressed sample, 
were dominating in this experiment. 

The definitive confirmation of validity of formula (4) 
appeared only recently, when the paper [1] was pub- 
lished. Its authors admit the discrepancy between their 
model of magnetic anisotropy arising at the interface be- 
tween two magnetic layers and the experimental data. 
Specifically, they expected a quadratic sinusoidal angle 
dependence but observed additional peaks at <p — 90° 
and if = 270°, see Fig. 1 in this paper. Their model 
mimics quite well the major part of data and curves pre- 
sented in Fig. 4c [1] but fails to explain the presence of 
those mysterious 'additional peaks.' 

4. Discussion 

Rippled surfaces are well known in experimental 
practice. Some studies were already performed aim- 
ing to gain the full control on ripple formation pro- 
cess on various substrates [4] (sapphire), [5, 6] (sili- 
con), [7] (ZnO), or to investigate the influence of ripples 
on various physical properties, most notably the mag- 
netic anisotropy, exchange bias [8], or morphology of 
magnetic domains. Recently many papers are devoted 
to rippled surfaces of diluted magnetic semiconductors 
(DMS), with (Ga, Mn) As being probably the most fre- 
quently studied substance in this class [9, 10]. The ac- 
tive area of research, both theoretical [11,12] and exper- 
imental [13] are competing anisotropics, uniaxial and 
tetragonal, present in thin layers of this compound. 

The rippled surfaces were approximated in literature 
in many ways, usually as a train of sinusoidal waves, 
or as a periodic series of Gaussian-shaped peaks or as 
a periodic set of flat islands. Here we propose yet an- 
other approach, namely the rippled surface may be seen 
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as being build by many identical, infinitely long half- 
cylinders, aligned parallel to each other. Obviously, the 
period of such a structure is equal to 4p, where p is, 
as previously, the individual cylinder's radius. Addition- 
ally, we neglect eventual interactions between cylinders. 

We test our theory, given in sec. 2. Using scanned 
data from Ref. [1], we try to fit them to the expression 

E{0) = a + K s \cosO\ + K v ,cos 2 O, (5) 

i.e. taking into account only the surface anisotropy and 
conventional uniaxial anisotropy. The constant a is ir- 
relevant, but has to be fitted in order to simulate ex- 
perimental data correctly. It is therefore not reported in 
Tab. 1, where the results for two available samples, with 
different thicknesses f, are shown. The ripple's period 
was reported in [1] as being equal to 22 nm, so the es- 
timated mean radius of curvature is 5.5 nm. This should 
be compared with the diameter of microwires used in 
[2]. Looking at the relation (3), it is easy to see why the 
surface anisotropy term could not be detected in ear- 
lier experiments: now the squared radius of curvature is 
some 3.0 x 10 5 — 2.25 x 10 6 times lower, and so many 
times the expected magnitude of the effect should in- 
crease. The reported uncertainties for both anisotropy 



Table 1 : Fitted parameters of expression (5) 
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constants, K s and K u , are most likely seriously underes- 
timated, some 2 — 3 times, by our quick and dirty fit. It 
is quick and dirty because the fitting procedure has no 
information concerning the uncertainties of individual 
measurements, and, consequently, treats all the data as 
being exact. Even the discretization errors, being a re- 
sult of manual scanning procedure, go unattended. 

Despite of these deficiencies, the trend is clear: the 
magnitudes of both anisotropy constants slightly de- 
crease with increasing sample's thickness. The decrease 
of K s probably has its roots in decreasing height of rip- 
ples, while their period stays unchanged during sam- 
ple growth, hence the curvature radius p effectively in- 
creases. This is probably also the reason for evidently 
'rounded' shape of peaks visible in Fig. 2. This feature 
may be also explained by finite lengths of individual 
cylinders, their misalignment, or even weak, but long- 
range, interactions between them. Since it not visible so 
clearly in Fig. 1, then we should attribute it to broad- 
ening of height (and, consequently, of curvature radii) 
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Figure 1 : MAE for the sample 2.0 nm thick. In addition to the best fit- 
ted line (blue) shown is the 'pure' uniaxial part (scaled) of anisotropy. 
Note the remarkable difference between the two near <p = 90° and 
ip = 270°. 156 measurements. 



distribution while the sample gets thicker. On the other 
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Figure 2: Same as Fig. 1, but for sample 5.2 nm thick. 297 data points. 
For clarity, the 'classical' uniaxial anisotropy term is not shown. 

hand, the drop in value of ordinary uniaxial anisotropy 
K u originates most likely from the strain relaxation far 
from possibly mismatched substrate. 

It is doubtful whether the presence of sharp, non- 
differentiable features, existing in reality on any ex- 
perimental curve will ever be possible to convincingly 
demonstrate using the data alone. O'Grady [14] shows, 
using symmetry arguments, that angular variation of 
many magnetic properties may be described either as 
even or as odd series of cos 8. Similar ideas, related to 
variation of coercivity or exchange bias field, were pre- 
sented even earlier [15]. Here we show one more pos- 
sibility: the surface anisotropy is described by a sin- 
gle | cos 1 term, rather hard to approximate by only few 
terms of even cosine series. 
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It remains to be explained why in Tab. 1 is ex- 
pressed in kJ/m 3 rather than in kJ/m 2 . This is intended, 
as it illustrates convincingly (in last column) compa- 
rable shares of both types of anisotropy in free en- 
ergy (not its density!). In fact, what we present there 
is the quantity K s = K' s /t, where K[ is the true sur- 
face anisotropy, expressed in kJ/m 2 , as it should, and 
t is the sample thickness. Taking this into account, we 
have K' S (2.Q nm) = -0.36 erg/cm 2 and K' s (5.2 nm) = 
-0.84 erg /cm 2 , respectively. One may wonder why the 
two estimates differ so much. It is less surprising when 
we compare the samples' thickness (2.0 and 5.2 nm) — 
in both cases smaller than the radii (5.5 nm) of our hy- 
pothetical half-cylinders. This means that our cylinders 
must be far from perfect, they are most likely flattened, 
what certainly affects their curvature radii, p. Neverthe- 
less, the values of K s decrease with sample thickness, as 
expected. Published values of \K' S \ are scarce, ranging 
from 0.032 erg/cm 2 to as high as 1.17 erg/cm 2 for Fe 
deposited on GaAs [16]. Our result is of the same order 
of magnitude. 

Let us now estimate the exchange energy per single 
Fe-Fe pair. The density of elementary cells on (001) 
surface of a-Fe is n = 1/a 2 « 1.214 x 10 19 rrT 2 , where 
a = 2.870 A is a-Fe lattice constant. Therefore the ex- 
change energy, E ex , per a x a square element of a sur- 
face is K' s /n, i.e. -2.965 x 10~ 23 J for thinner, and 
-6.853 x 10~ 23 J for thicker sample. From formula (3) 

we get J = £ex(^f, that is -3.937 x 10~ 21 J and 
-9.099 x 10~ 21 J, respectively, when S = 2.22 [^ B ] 
and S = a V3/2 * 0.215 nm. Those values should still 
be divided by the number of nn Fe pairs (4) residing in 
a said a x a surface element. This is because the nearest 
neighbor for any given Fe surface atom is the one lay- 
ing deeper, inside the elementary cell - as pure iron has 
bcc structure. This fact has been already accounted for 
by expressing 5 {nn spacing) as an appropriate fraction 
of the lattice constant. Finally we obtain 7(2.0 nm) = 
-0.98 x 10~ 21 J, and 7(5.2 nm) = -2.28 x 10~ 21 J. For 
comparison, Ref. [17] quotes J = -1.21 x 10~ 21 J 
for pure a-iron. The correspondence is amazingly good, 
especially that our model completely neglects RKKY- 
type exchange, certainly present there, and the estimate 
is made as if the surface was perfectly flat. 

5. Conclusions 

The surface anisotropy form, presented here, seems 
to explain the observed features of magnetic anisotropy 
energy simply formidably. The shape of angular depen- 
dence of MAE is reproduced much better than by any 



other model. The deduced values of nn exchange cou- 
pling strength are in good agreement with those ob- 
tained independently. Moreover, they are in full accor- 
dance with intuitive understanding what makes the sur- 
face layer: no more than two crystal planes are involved. 
Consequently, the surface layer thickness is lower than 
the size of a unit cell. Yet, such effect can be eas- 
ily observed only at nanoscale, that is in samples thin 
enough. Only then its magnitude is comparable with 
ordinary uniaxial anisotropy (see the last column of 
Tab. 1). One may expect that, at least in the case of iron, 
a ~ 500 nm layer is thick enough to effectively mask 
surface anisotropy effects. 

It is amazing that our original, idealized model of 
non-interacting, infinitely long half-cylinders, works so 
well. It is likely that the presence of elongated, but fi- 
nite length structures, present on nominally flat sur- 
faces, even those obtained by MBE technique, is suf- 
ficient to generate this form of anisotropy. On the other 
hand, it is doubtful whether it will ever be used to deter- 
mine some parameters it depends on. It is because the 
presented surface anisotropy term is rather sensitive to 
the fine details of a surface. Those are probably easier to 
investigate using one of microscopic techniques. Never- 
theless, using its peculiar angular behavior, and treating 
it as a 'background' of known shape, one should be able 
to determine important material's parameters with bet- 
ter accuracy than it was possible earlier. 

The presence of non-negligible surface anisotropy, 
generated by surface curvature, in addition to the edge- 
related effects, will affect the operation of future spin- 
tronic devices. 
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